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Diffusion properties of a self-avoiding polymer embedded 
in regularly distributed obstacles with spacing a = 20 and 
confined in two dimensions is studied numerically using the 
extended bond fluctuation method which we have developed 
recently. We have observed for the flrst time to our knowl- 
edge, that the mean square displacement of a center monomer 
(t>M/2{t) exhibits four dynamical regimes, i.e., 4'M/2it) ~ f^"^ 
with Vm ~ 0.6, 3/8, 3/4, and 1 from the shortest to longest 
time regimes. The exponents in the second and third regimes 
are well described by segmental diffusion in the "self-avoiding 
tube". In the fourth (free diffusion) regime, we have numer- 
ically confirmed the relation between the reptation time rj 
and the number of segments M, r^ oc M^. 



I. INTRODUCTION 

The reptation theory, proposed for concentrated so- 
lution of linear polymers |l[, is based on the assump- 
tions that a polymer moves in a "tube" made of the 
other polymers and that the tube does not deform till 
the polymer of interest creeps out of it. Since the the- 
ory was proposed, various investigations have been done 
on dynamic properties of a polymer in networks of im- 
mobile obstacles. As for the case of concentrated solu- 
tion, Baumgartner et al y] studied dynamics of a freely 
jointed chain with the Lennard- Jones potential in a dense 
network of frozen-in chains. Evans and Edwards ||^, on 
the other hand, proposed an efficient lattice model with 
stochastic local jump motion of segments in regularly 
distributed obstacles on the cubic lattice, which simu- 
lates such phenomena as dynamics of a DNA molecule in 
agarose gel. 

Study of diffusion property of a single real polymer 
(self-avoiding chain) in fixed and low density network of 
obstacles is of importance, since details of DNA dynam- 
ics in agarose gel, etc. ^-^ have been studied by many 
experimental works along with recent developments in 
technical aspects |7|§. To see entanglement effect, in- 
vestigations of a self-avoiding polymer in two dimensions 
would be desired. In two dimensions, a molecule with 
relatively low polymerizations can entangle with many 
obstacles since it extends larger than a polymer in three 
dimensions. 



Recently, Maier and Rodler observed conformation 
and diffusion of a single DNA molecule electrostatically 
bound to fluid cationic lipid bilayers by fluorescence mi- 
croscopy. They measured radius of gyration Ri and ob- 
tained the relation {B?j) ~ N'^'^ where N is the number 
of base pairs of the molecule and v — 0.79 ± 0.04. For 
diffusion coefficient of the center of mass Dq they showed 
Dg oc A'^"^. These results agree well with the theoretical 
predictions for a self-avoiding polymer in two dimensions 
without hydrodynamic interactions. It seems feasible to 
arrange obstacles on that layer and somehow measure 
displacement of a particular part of DNA to examine 
diffusion property of a self-avoiding polymer in a fixed 
network of obstacles in two dimensions. 

In the present work, dynamic properties of a single 
polymer in regularly distributed obstacles in two dimen- 
sions is studied by means of the extended bond fluc- 
tuation method (e-BFM) which we have developed re- 
cently jl^]. From the simulated mean square displace- 
ment of a center monomer 4'M/2{i)i we can distinguish 
clearly four dynamical regimes, as is the case of a Gaus- 
sian chain ifllj in the same network of obstacles. The 
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3/8, 3/4, and 1 from the shortest to longest time regimes. 
The exponents in the early three regimes differ from those 
of the Gaussian chain, but they are well described when 
the excluded volume effect is taken into account. In par- 
ticular, dynamics in the second and third regimes is inter- 
preted as segmental diffusion in the "self-avoiding tube," 
where displacement As along that tube corresponds to 
aAs"^ in real space with a and v being the tube diameter 
and the exponent relating radius of inertia Rj and the 
number of segments M as Ri ^ Af, respectively. The 
reptation time t^, as for the crossover time between the 
third and fourth (free diffusion) regimes, is ascertained 
to be proportional to M^ with M being the number of 
segments of the chain. 

In the next section we explain the details of the simula- 
tional methodology and parameters which we have used. 
The scaling theory of tube model for a self-avoiding poly- 
mer in the network of fixed and low density obstacles is 
described in the Sec. [II. We then present the results and 
their scaling analyses compared to the theory in Sec. IV, 
and conclude them in the final section. 



II. METHOD 

The bond fluctuation method (BFM) ^j^ is known 
to provide effective Monte Carlo (MC) algorithm to study 
static and dynamic properties of various phenomena in- 
volved in polymers |1J]. We have adopted a square lat- 
tice version of the method where a polymer is described 
by a sequence of monomers which occupies four sites of 
the smallest square cell. The monomers are not allowed 
to occupy the same site (excluded volume interaction or 
effect (EVE)) and connected by bonds which can vary 
its length I in the range 2 < I < \/T3. The elemental 
MC step is a trial of displacement of each monomer by 
one lattice spacing to a randomly chosen direction. The 
process satisfies the self-avoiding condition. 

The extended-BFM which we have proposed pOJ and 
used in the present study incorporates a non-local move- 
ment of the 's-monomer' into the BFM. The s- monomer is 
defined as the one whose nearest neighbors separate not 
greater than vl3 with each other. An s-monomer can 
be displaced to any monomer pairs between the nearest 
s-monomers by a stochastic process which fulfills detailed 
balance condition. This new non-local process can be in- 
troduced in any relative frequency into the conventional- 
BFM above mentioned. We have used the following com- 
bined process: at even (odd) time step of the c-BFM, 
the new process for even (odd) numbered s-monomers is 
tried. 

The e-BFM has been initially introduced to overcome 
the trapping or pinning difficulty which the c-BFM is 
faced in gel electrophoresis under a large electric field. 
It has turned out that the method reproduces static and 
dynamic properties also under a vanishing field correctly, 
and more faster than the c-BFM for polymer dynamics 
in a space with fixed obstacles. In fact, the ratio A of 
the diffusion constants obtained by the c-BFM and e- 
BFM is less than unity: A ~ 1/2 and 1/3 for a = oo (no 
obstacles) and 20, respectively, with a being the distance 
between the obstacles. 

By means of the e-BFM we here mainly observed 
the mean square displacement of the center monomer 
4>M/2{t) a-s a representative of inner monomers. The ini- 
tial configurations are provided by the self-avoiding walk 
algorithm where monomer sequence is made by iteration 
of trials of making a new end. If it finds bond cross- 
ing or monomer overlap, it stops making the chain and 
go to another starting position. Although this algorithm 
can give independent equilibrium configurations, the av- 
eraged time needed for a complete configuration turned 
out to be quite large due to the low acceptance rate for 
each trial. So we have relaxed the condition for accepting 
a new end as that it allows up to uj trials. Here we have 
chosen w = 7 ~ 14. The number of samples given at 
t = 0, the start time of the idle run, was 512 ~ 788. The 
interval of the measuring run tf — U was taken so as to 

[4'M/2itf - U)] ^^^ ^ CRi, with C > 2. The system size 
was 3Af X 3M in which even a fully extended polymer 



cannot interact with itself through the periodic bound- 
ary conditions introduced. We distributed, regularly in 
a spacing a = 20, obstacles of a unit square cell, with 
which any monomer site cannot overlap (the EVE). 

Let us here compare some results simulated by the c- 
BFM and the e-BFM. As an example, we show (/)M/2{t) 
for a = 20 and M = 150 in Fig. |l|. At t < Te, t^ being 
the crossover time between regimes I and II introduced 
in the next section, (j){t/Aoc) with A^a = 0.5 simulated 
by the c-BFM and (j)M/2{t) by the e-BFM coincide with 
each other, and exhibit the power-law growth which is 
proportional to i°'^. The latter is the characteristics of 
monomer diffusion in free space. At t > Te, on the other 
hand, where monomer diffusion is affected by the sur- 
rounding obstacles, <j>M/2{t) by the e-BFM is seen to co- 
incide with (/)jvf/2(i/^2o) with A20 = 0.3 by the c-BFM. 
The factors Aoo and A20 are just the ratios A of the diffu- 
sion constants mentioned above. Thus the results shown 
in Fig. [l| demonstrates that, there is no fundamental dif- 
ference in diffusion properties of the polymer simulated 
by the c-BFM and the e-BFM. 
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FIG. 1. The mean square displacements of the center 
monomer of the M = 150 chain in regular distributed (spacing 
a — 20) obstacles. The symbols D, Q, and A are of the 
e-BFM, c-BFM for t <Te, and c-BFM for t > Te respectively. 



III. SCALING THEORY 

We begin with brief summary of the scaling argument 
on the tube theory of a Gaussian chain moving in a 
network of fixed obstacles ||ll|. There exist four time 
regimes in which the chain exhibits different dynamics. 
In the earliest regime denoted by regime I, the n-th in- 
ner monomer ffuctuates without perceiving wall of the 
tube, and its mean square displacement (f>nit) is writ- 
ten as (f>nit) ~ t'^'" with Vm = 1/2 |l£|. Here we put 
C/ksT ~ 1 with C being the friction coefficient of a seg- 
ment and the segment length b unity. In regimes II and 
III, the polymer motion perpendicular to the tube coordi- 
nate is restricted. Diffusional displacement of the prim- 



itive chain along the tube (As^) in a time interval t is 
written as (As^) ~ at^'"^ and at^ before (regime II) and 
after (regime III) the Rouse relaxation time t^ ^ A'P, 
where a is unit length of the primitive chain which is set 
equal to the tube diameter. Such displacement along the 
tube corresponds to that in real space (Ax^) through 



shown schematically in Fig. 0. Each blob of an average 



(Ax^) ^ (As^ 



(1) 



with v = 1/2 since the tube is Gaussian. Thus we ob- 
tain Vm = 1/4, 1/2 for regimes II and III, respectively. 
In regime IV, the whole chain creeps off the tube over 
the length scale of Ri , and motion of the segment coin- 
cides with that of the center of mass which follows the 
linear relation </>„(i) ~ Rjt/Td- Here t^ ^ M^ /a? is the 
reptation time. 

In regime IV of a model simulating a Gaussian chain in 
the regularly distributed obstacles [^ , viscosity yy, which 
is compared to r^ |ll|, was measured up to L = 50 
with L being the number of beads and was shown to 
be proportional to 7^3.4i±o.i4 ^gj^g 7^ > 20 jl^. Recent 
study on this model of the most concentrated case clearly 
showed Vm = 1/4 in the shorter time region, i.e. regime 
II (regime I was not seen) [^ . 

Now let us turn to dynamics of a real chain (with the 
EVE) in space with fixed obstacles. In regime I, as for 
a Gaussian chain, inner monomers fluctuate without the 
tube constraint but with the EVE. In this case, or with- 
out obstacles (a = 00), it is known that (j)M/2{t) oc t^ at 
t < Tr, where z = 1/(1 + {1/2vf)) with up = 3/{d + 2) 
and d being dimension of the space, and r^ is the rota- 
tional relaxation time for a self-avoiding chain pq] . 




FIG. 2. The schematic representation of a chain of blobs. 
The polymer configuration drawn is one of instantaneous ones 
observed in our simulation with M — 400. The blobs are 
drawn intuitively from the polymer configurations averaged 
over time interval of about Te. In an instantaneous configu- 
ration there some local configuration of hernia type are seen. 

A polymer configuration in regimes II and III is con- 
sidered as a chain of blobs along the tube coordinate as 



size Ob consists oi g {- 



Ijvp 



) monomers on average. In 



terms of the tube coordinate, the end-to-end distance R 
is written R/ay, ^ (M/g)"'', which impHes Ri ex M"'' for 
the original chain. In fact the latter has been confirmed 
numerically at least for low concentration of obstacles in 
two dimensions as shown in Fig. 0. Thus the EVE plays 
an essential role in determining a "self-avoiding tube" in 
these regimes. 
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FIG. 3. The relation between radius of gyration Ri and 
length M at a = 20. 

For dynamical aspects under a given polymer config- 
uration (and so a self-avoiding tube), however, the EVE 
between the blobs of length scale Ob is considered to be 
screened out |19|-pif|. Accordingly, the diffusion process 
of blobs along the tube is the same as that of the Gaussian 
chain [n|. Thus, the exponents of the mean square dis- 
placement of the center monomer along the tube (t>M/2\\ (t) 
are given by 1/2 and 1 for t < tq and t > Tg, respectively, 
where Tg — Te{M/g)'^ is the Rouse relaxation time of the 
chain of blobs. Now taking into account eq. (0) with 
V = vp, we obtain the exponents Vm = j^f/2 and vp in 
regimes II and III, respectively. In regime IV Vm = 1 as 
in the Gaussian chain. 

Now we derive <j)n{i) with n ~ M/2 in each regime 
more quantitatively. As the Rouse and the reptation the- 
ory claim relations Tg/Te ^ [M/g] and Td/Te ^ [M/g] , 
the pre- factors of 4>n{t) in the all regimes can be deter- 
mined by matching its scaling forms in the neighboring 
regimes of Tg, Tg and r^. Setting b, (/ksT = 1 as before 
and noting that 0„ ~ R^ [t/'''d] in regime IV, we obtain 
0„ (t) and the cross-over times as listed on Table H. 
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TABLE I. Expressions of (j)n{t) in the four regimes 



Lastly we consider the relaxation time td of the auto- 
correlation function of the end-to-end vector defined by 
P{t) = (R(t)R(O)) oc exp[-t/T£,]. For regime IV the 
reptation theory [1^ predicts that tjj, called the rep- 
tation or disengagement time, is obtained from one- 
dimensional diffusion equation of the primitive chain. It 
is given by td — L^/ir'^Dc, where Dc — l/M is the diffu- 
sion constant along the tube, and L is the length of the 
primitive chain and is given by aL — M for a Gaussian 
chain ( note we put ksT/C, — 1 and 6 = 1). We can 
rewrite t^ for a self avoiding chain by replacing L by 
UhM/g. This is because the motion of the chain along 
the tube is described by dynamics of a sequence of blobs 
which behave as a Rouse chain. Therefore T]j and t^ in 
Table | are related as 
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(2) 



for a real chain. 



IV. RESULTS 

In Fig. ^ we show the result of the mean square dis- 
placement of the center monomer 4>Ml2{t) obtained from 
our longest chain M = 400 in two dimensions. As listed 
on Table 1 in the previous section, (t>M/2{t) is expected 
to exhibit different power-law dependences on t in the 
four regimes: with d = 2 and so vp = 3/4 and z = 3/5, 
i^m = 3/5,3/8,3/4 and 1 in from regime I to IV, re- 
spectively. These theoretical results are well reproduced 
in our simulated results in the figure. The auxiliary 
line in the shortest time interval (regime I) stands for 
the least square fit of our separately simulated data of 
the M — 400 chain with a — oo (no obstacles) to the 
power-law function of t. The exponent is estimated as 
Vm = 0.60 ± 0.01. The data of a = 20 lie upon this hue 
at 10 < i < 2 X 10^. as well. 



10* 
10^ 
10^ 

10^ 

IQi 

10° 

10° 10' 10^ 10' lO'* 10^ lO*" 10^ 10*^ 10'^ 

t 



7 


a- 
M 


, ,, , , ,, , . ,, . , 
:20 

=400 


~y 


/ 






3 ^ 




~i?? 


7 




^-"^ 




- 


; ~o.^ 


y 




1 


'e 


_ 




. .1 . 


. .1 . . .1 . . .1 . . 




1 . . . 



FIG. 4. The mean square displacement of the center 
monomer of the M = 400 chain in regular distributed (spac- 
ing a — 20) obstacles. 

The tube model predicts the crossover time between 
regimes I and II, where fluctuating monomers start to 
feel the wall of tube and to form blobs, as (i, 4>M/2{t)) = 
(te, a^) . From inspection of the figure, its ordinate seems 
significantly smaller as compared to a^. It is estimated 
as Ob ~ a/2. 

We may think of another diameter dx, which corre- 
sponds to diameter of a straight tube with hard wall. 
The mean square displacement of a monomer perpendic- 
ular to tube 4>M/2±it) for i ^ Tg is estimated as 



H//2±(i> Te) ~ 



dT-2 dr-2 

dT-2dT-2 

EE 

1=1 3=1 

1 [dr - 1] [dr - 3] 



(3) 



Replacing 4>M/2±it) by 0M/2(Te) = < = 114 ± 5 and 
substituting it to eq. (|), we get dr ~ 28.2 ± 0.57. This 
indicates that the tube diameter dx is bigger than a since 
the "tube" is not a tube with hard wall but is formed in 
a space of the point obstacles. 

The data in Fig. ^ also demonstrate nicely the chain 
dynamics in the "self-avoiding tube" described in the pre- 
vious section. In fact the predicted exponents 3/8 and 
3/4 in regimes II and III respectively are ascertained 
with the appropriate crossover point {t,(j>M/2{t)) — 
{Tg,a\)Rj). Finally at i > r^ crossover from region III 
to IV occurs, and the linear relation 4>M/2(t) ^ i is re- 
produced in regime IV. 

In order to confirm further the compatibility of theo- 
retical arguments on the polymer dynamics in the self- 
avoiding tube, we have carried out scaling analyses on 
4>M/2{t) with different M {— 100 ^ 400). Some raw data 
are shown in Fig.pl 
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FIG. 5. The simulated data of (f)M/2{'t) with M = 100, 
200, 300 and 400 at a = 20. As expected from the fact 
Ri > ab for these M they coincide with each others in regime 
I and a shorter time range of regime II. 

In Fig. H double logarithmic plots of 4>M/2it/Te)/ahRi 
in terms of t/rg around crossover between regimes II and 
III are shown for M = 100, 200, 300 and 400. Here 
the value tq is determined as Tg — t^ [Ri/a\^] , using 
Rj and fp in Fig. |3| and Tg and at, from Fig. ^. It can 
been seen that the scaling fit is satisfactory for M > 
200. The auxiliary lines fitted to the data yield exponents 
0.38 ± 0.01 and 0.73 ± 0.02 in the regimes II and III, 
respectively, and they cross at the point (1,1) within our 
numerical accuracy. The data of the shortest M = 100 
chain, on the other hand, deviates from the others in 
early time region t/rg ^ 2 x 10^^, where it is in the 
regime I. So it can be said that dynamics in regimes II 
and III for M > 200 is interpreted as monomer motion 
through blobs along the self-avoiding tube. 



10^ 



J3 



^ 



10° 



10" 



M=100 B 
200 - - 
300 
400 




10" 



lO'^ 



10' 



t/Te 



FIG. 6. The double logarithmic plot of cj)n[/2{t/fe)/ahRi 
with t/Tg. 



We show double logarithmic plots of 4>M/2{t/'i'd)/ R^i 
against t/r^ in regimes III and IV in Fig. 0. Here r^ = 

Te [Ri/a\^] '"'^ (ex A'P). The data in each regime are 
scaled to a single curve for M > 200, and they cross 
near (1,1) as well except for the shortest (M = 100) 
chain. This result supports the crossover scenario from 
the reptation regime of HI to the over-all diffusion regime 
of IV. 
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The double logarithmic plot of 4>M/2{'t/Td)/ Ri 



Next we discuss behavior of some other quantities than 
4'M/2{t) niainly in regime IV. Figure B shows the semi- 
log plot of P{t) for M = 50 ^ 400. At long time regions, 
which corresponds to regime IV with 0m/2(O '^ ^i Pi^) 
of each M is well fitted to an exponential function. The 
plots of the relaxation time extracted from P(t) and the 
one evaluated by means of the r.h.s. of eq. (0) against 
M are shown in Fig. g. The two estimates of td co- 
incide with each other even up to numerical coefficient 
and clearly display the slope 3 within error of 3%, i.e. 
3.0 ± 0.1, for M > 70. The result strongly supports the 
argument that the primitive chain, whose unit length is 
identical to blob size Ob, diffuse along itself. 
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FIG. 8. The autocorrelation function P{t) of the 

end-to-end vector observed for each lengths. 
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FIG. 9. The rotational relaxation times td obtained from 
Fig.|andJ^[^]''-/V,ofeq. §. 

In such a free diffusion regime monomer motion is iden- 
tical to that of the center of mass. The pre- factor of linear 
fit of the mean square displacement of the center of mass 
'/'c.m.(i) in the longest time range (which includes regime 
IV) gives us diffusion constant Dq mU. Typical data of 
(i) are shovirn in Fig. 10. We have determined Dq 



by the linear fit of the data in t ^ Td, and plotted them 
against M in the inset of Fig. |lO|. 

The exponent a for Dg ^ M" is estimated as a ~ 
— 1.5 for large A/. It is different from what Gaussian 
tube theory predicts (a ^ —2). It is, instead, readily un- 
derstood if we plot Td ( oc M'^ ) and RJ/Dg(oc AP"-" 
) against M simultaneously. As seen in Fig. |l^ the scal- 
ing relation (3 — 2v ~ ol holds as is derived from scaling 
argument of the reptation theory |2^ . 

We have obtained a ~ —1.5 using the data with 
M = 70 - 400. When the data with A/ = 30 - 150 
are used, on the other hand, it becomes a ~ —1.7. This 
is consistent with the results obtained in our previous 
work [nOl and is quantitatively in agreement with other 
results such as reported in [|l6[. Such a large magnitude 
of a obtained from smaller M may be attributed to tube- 
length fluctuation which gives rise to a leading correction 



to the scaling described so far. Actually we have con- 
firmed /3 = 3.2 ± 0.1 from the data with M = 30 - 150 
as seen in Fig. g. Then the scaling relation ^ = 2v — a 
still holds for these exponents as shown in Fig. ^ 

Now we taking into account such correction, the devia- 
tion of the M = 100 data from the scaling curve in Fig. 
is naturally interpreted that M = 100 is in the crossover 
range of the parameter space of M . 
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FIG. 10. The mean square displacements of the center 
of mass cf>c.-m.{i) and their diffusion constants Dq, calculated 
from tangent of final linear regime, are plotted against M. 
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FIG. 11. The scaling relation f) — 2v — a. 



V. CONCLUSIONS 

We investigated on the dynamics of single long self- 
avoiding polymer in the regularly distributed [spacing 
a — 20, i.e. concentration c = 1%] obstacles in 2D by 
the e-BFM and analyzed the results especially in con- 
nection with the tube model. The following four regimes 
are observed in the mean square displacement of cen- 
ter monomer 4>M/2{t) c< t^'^ as v^ ~ 0.6, 3/8, 3/4 and 



1. These exponents are explained by the "self-avoiding [21] P.-G. de Gennes, Scaling Concepts in polymer physics 

tube" model. That is, the defect dynamics through blobs (Cornell University Press, Ithaca, New York, 1979). 

along such a swelled tube give the exponent of I'm as 3/8 

and 3/4 at the second and the third regimes respectively. 

We have shown the validity of such argument through 

the scaling analyses on the crossover regions and the free 

diffusion regime. 
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